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SUMMARY 


c The problem of the build-up of lift on two— and three-dimensional 
wings flying at high speeds is discussed as a boundary—’ value problem for 
the classical wave equation. Kirchhoff*s formula is applied so that the 
analysis is reduced, just as in the steady state, to an investigation of 
sources and doublets. Some simple applications of this method are con- 
sidered, including the determination of the starting lift of a three- 
dimensional wing and the potential functions for some types of unsteady 
vortex motion. 


INTRODUCTION 


The usual idealizations introduced in the development of linearized 
aerodynamic theory describe a frictionless, perfectly elastic, model 
fluid. As is well known, the effect of small disturbances in such a 
fluid can be analyzed by means of the familiar wave equation which, in 
terms of the perturbation velocity potential cp, can be written 

'Pxx + ^yy + ^zz = a 2 

o 

where 

a Q speed of pressure propagation 

t * time 

x,y,z spatial coordinates 

Equation (l) applies to a flow field which is stationary at large dis- 
tances from the disturbance region; furthermore, the axial system is 
stationary relative to the fluid infinitely distant so that, if a moving 
wing is being analyzed, the wing moves with respect to the x,y,z axes. 


'Vt' 


( 1 ) 
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It is possible to use a transformation which will bring the wing to 
rest with respect to a new axial system, but this process leads generally 
to a different equation for the perturbation potential. The case of a 
wing moving at subsonic speed constitutes an exception, for here the 
Lorentz transformation, familiar in modem physics, will fix the origin 
of the axial system on the wing while leaving the wave equation invariant. 
Such is not the case when the wing travels at supersonic speeds, so that, 
if it is desired to fix the axial system to the wing, the wave equation 
must be abandoned. This latter alternative has been the more extensively 
studied in aerodynamics, but, because of the large body of knowledge con- 
cerning equation (l) which is available in mathematical physics, the 
present report is based upon the wave-equation approach to problems of 
unsteady motion. 

The solution to equation (l) known as Kirchhoff’s formula (see ref- 
erence 1) is found to be of considerable use in unsteady-motion problems 
involving thin wings with supersonic edges. The problem is reduced to 
one of summing elementary solutions, analogous to sources and doublets 
in steady flow, over a region determined by the position of the wing as 
well as it 3 traversed path. The theoretical development leads naturally 
to the concepts (defined later) of, first, inverse sound waves, which 
have a counterpart in the Mach forecones used in steady-state wing theory; 
second, acoustic plan forms; and, third, homogeneous flow, which reduces 
in part to the familiar conical flow as the wing approaches a steady 
supersonic velocity. 

In connection with boundary— value problems involving noninteracting 
surfaces, a theorem will be given which makes it possible to build up a 
three-dimensional supersonic— edged wing by superposition of elements of 
two-dimensional wings. This result is an extension to unsteady flow of 
a theorem given previously by Lagerstrom and Van Dyke (reference 2). By 
means of this theorem, problems involving super sonic-edged wings per- 
forming any prescribed maneuver can be solved, provided only that the 
behavior of a two-dimensional wing executing the same maneuver is known. 
Thus, the relatively large amount of material available for two-dimensional 
unsteady motion at supersonic speeds can be carried over directly to the 
supersonic— edged wing. 

In the last part of the report, some simple applications of the gen- 
eral methods are considered, including the determination of the starting 
lift of a three-dimensional wing and the potential functions for some 
types of unsteady vortex motion. 


UST OF IMPORTANT SYMBOLS 


a Q free— stream speed of sound 

c 0 wing root chord, or maximum chord of a wing 
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c 


c la 

C L 


L 

Mo 

q 

q 

t* 

t 

v 0 

w 


x,y,z 


a 

P 

X 

p o 

<P 

A 


u 


wing chord 

section lift-curve slope 


lift coefficient 

coefficient of rolling moment due to yawing 


[ 


a Cj 

a(rt/2V 0 ) _ 


lift 

free— stream Mach number 

loading coefficient (pressure on the lower surface minus pressure 
on the upper surface divided by dynamic pressure) 

free-stream dynamic pressure ( I c ° v ° 2 ) 

time 
a 0 t ’ 

free-stream velocity 


vertical component of perturbation velocity 
Cartesian coordinates 
angle of attack 


J Jl -Mo 2 | 

slope of stream surface in free-stream direction 
free-stream density 


(*) 


perturbation velocity potential 

discontinuity in the quantity in question across the z = 0 plane 


Subscript 


value on upper surface (z = 0 plane) of a wing 
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THEORETICAL DEVELOPMENTS 


Basic Equation and Solutions 


The partial differential equation which forms the basis of the pres- 
ent study has been written as equation (l). It is possible to put this 
equation in a more convenient form by introducing the notation 

t = a Q t» ( 2 ) 


bo that the dimension of t is length just as are the dimensions of the 
geometric variables x, y, and z. Equation ( 2 ) y together with equation 
(l), yields the canonical form of the wave equation 


<Pxx + Oyy + ‘Pzz = * tt 

and it is this form which will be considered. 

The first task is to study the relation between the motion of the 
wing and the coordinate system. As has already been mentioned, equa- 
tion ( 3 ) is valid for a flow field produced by a wing moving relative to 
a fixed coordinate system. It is pertinent to consider the possibility 
of finding a transformation which will (l) fix the origin of the axial 
system on the wing and, at the same time, ( 2 ) retain the wave equation 
as the governing equation of the flow. Certainly the first of these 
requirements is simple to fulfill if the second is neglected. However, 
because of the simplicity of the wave equation and, what is more impor- 
tant, because of the great amount of developmental study that has been 
expended on it, the second requirement is not without justification. 

The following transformation (known in relativity theory as a Lorentz 
transformation, or as a hyperbolic rotation) 

x-Mpt 

7.1-Mq 2 

y 

Z 

t-MpX 

Vl-Mp 2 



t = 

n = 
i = 



MCA TN 2256 


5 


where M Q = — <1, will satisfy both conditions (l) and (2). For example, 

suppose that the wing is moving along the x axi3 with velocity V 0 . 
Application of . equation (4) to equation (3) makes the origin of the new 
axial system also travel along the x axis with a velocity V 0 . This is 
seen to be so since | is always zero when x = V 0 t* = Mot. Hence, the 
5 axis is "fixed" on the wing. As to the second condition, a straight- 
forward exercise in partial differentiation yields 

<P|| + ^qq + c ^ = 'Ptt (5) 

so that, in going from x,y,z to |,q,£ space, the wave equation 
remains invariant; consequently, both the requirements mentioned have 
been fulfilled. 

It is instructive to consider briefly the consequences of applying 
the Lorentz transformation. Although the wave equation remains invariant, 
such physical quantities as length and pressure do not. For example, a 
wing with a chord c 0 in the x,y,z space has, according to equation (4), 

a chord Cq/Vi-Mq 2 in the £,q,C space. Furthermore, the loading 
coefficient which, on the basis of linearized theory, is given in the 
XjyjZjt apace by 

Ap = - ^ (6) 

q VqMo St 

becomes for the |, q, £, t space 

Ap 4 { 1 Sq? Sq)^ (7) 

V 0 Vl- Mq 2 S t S| / 


If the wing motion is steady and there are no transient effects, equa- 
tions (5) and (7) are independent of time and together with the resulting 
length transformations become 


<Pgg +<Pqq + 0 

Ap 

q ~ Wi-Mq 2 


( 8 ) 


£ = 


X 

yi-Mo 2 


n = y, C. » z 


J 
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These are immediately recognized to he, respectively, Laplace's equation 
for incompressible flow and the familiar Prandtl-Glauert compressibility 
corrections. 1 

The preceding discussion has an important qualification, however, in 
the fact that the velocity of the moving axial system cannot exceed the 
speed of sound. A glance at equation ( 1 +) serves to verify this statement 
since that equation shows M 0 must be less than 1 in order that | and 
t be real for real x and t. In fact, it has been shown that there is 
no transformation which will fix the moving axial system in a wing travel- 
ing at a uniform supersonic speed away from the original fixed axes and 
still keep the wave equation invariant. Therefore, for analyzing a wing 
in supersonic flight, it is necessary to abandon one of the two proposed 
requirements: either the axial system cannot be fixed in the wing, or 

the field equation must be modified. The latter of these two alternatives 
has been studied by several authors (e.g., references 3, 4 , 5, 6) but it 
is the former which will be considered in the present analysis. Further, 
since the axes cannot be made to travel as fast as the wing, they will 
not be made to move at all and equation (3) will be adopted throughout as 
the basic equation. 

Having decided upon the form of the partial differential equation, 
the boundary conditions must next be established. For any given time 
these conditions are similar to those studied in steady-state thin— 
airfoil problems; namely, either that the given slope of the wing surface 
is proportional to the vertical induced velocity (p z over the region 
occupied by the wing in the z = 0 plane, 2 or that the prescribed surface 
pressure is proportional to the timewise gradient in velocity poten- 

tial over the same region. The addition of time simply means that this 
region moves about in the z = 0 plane in conformity with the known 
direction and velocity of the wing. 

The additional condition is imposed that the induced velocities fall 
to zero on the surface which is formed by the envelopes of the spherical 
sound waves originating from the surface of the wing. 

The solution to equation (3), subject to the boundary conditions 
just mentioned, can be expressed by a formula which may be regarded as 

lTile equation |= x / J 1-Mq 2 would first read |= (x-I^t 0 )/y i-jyjg 2 
where t 0 is a constant representing the time required fpr the motion 
to reach its steady state. However, the x coordinate can always be 
translated to any fixed position without affecting any of the equations 
for potential, loading, etc. Such a translation is assumed to have 
been made in equation (8). 

^he z = 0 plane is assumed to be the "plane of the wing"; that is, if 
the angle of attack were zero and the wing had no thickness it would 
lie entirely in the z = 0 plane. 
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requiring either the evaluation of a double integral or the solution of 
a double-integral equation, depending upon whether a boundary-value 
problem of first or second kind is considered. This solution is known 
as Kirchhoff's formula. (See reference 1.) It may be written in a form 
convenient for aerodynamic applications as follows: 


1 p p[ 1 a _/ . . d <pUi,yi, 0 ,t-H?)'] 

V = -^JJ [^ a 5^r ,p(li - yi ’ zi ’ w o ) “ A ^ ? J 43 <9) 

b a 

where the A indicates the Jump (value on the upper surface minus value 

on the lover surface when applied to (p or of the function in 

oz 


passing through the z ± -0 plane, r = ^(x— x x ) +(y— yi) +(z-Zi) , 

r Q = 7(x- Xl ) 2 +(y-yx ) 2 +z 2 , and where the area of integration S a will he 

discussed in more detail later. 

The terms in the integrand of equation (9) can he shortened ty 
introducing the following notation: 


— A q>(xi,yi,z;L,t-r ) = ~ A [—■ 
r Q ozi r o LoziJ 


a [~ dr d <P(xi,yi,0,W) dr p ^ d f <P 1 

A dz! r ” L^Zi dr r J dz dr \_r J 

(10b) 

In this notation, the subscript r in equation (10a) means that r is 
to be held constant in the differentiation, and the prefix A obviates 
the necessity of indicating that the functions considered are to be 
evaluated for z x =0, since it indicates that the difference of the values 
of the function across the z x =0 plane is to be taken. The right-hand 
side of equation (10a) can be recognized as a term representing a source 
located in the Zi=0 plane, and the right-hand side of equation (10b) is 
seen to represent a doublet located in and with axis normal to the zx=0 
plane. The brackets [ ] about the functions in equations (10a) and 
(10b) have a special meaning which is defined in the following way: if 
f is a function whose value at a fixed point P depends upon the coor- 
dinates xi,yi,zi,t of a moving point Q, so that 

f = f(xi,yi,zi,t) 

then 

[f] = f(xi,yx,zi,t^) (11) 
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where r is the distance from P to Q. As an example, consider the 
potential <p at a point P due to a moving source, the location of which 
at any time is Q. Then cp satisfies the condition just mentioned that 
it depends on the coordinates x 1 ,y 1 ,z 1 ,t of Q. The brackets [ ] indi- 
cate that the potential [ cp ] depends not upon the source strength now at 
"time" t, but rather upon the source strength that existed "time" r 
ago. 3 For convenience, [cp] is referred to as the retarded value of cp. 


The expression for a doublet, equation (10b), is usually expanded 
as follows : 



I-® 

J 

dro 

» a 

.j i 

[*!?) 

dz dr 

l r J 

dz 

5 r oV r c 

)/ dz r o 

Sr J 


dz V r o 

Finally, equation (9) becomes 


J/i) [A<P] + 1 ar » 


9 = — r— 
4« 



1_ 

r 0 


dA<p 

Ldzi J 




*0 dz 


1_ 

r Q dz 




dAcp 

St. 


( 12 ) 


dS ( 13 ) 


The application of equation ( 13 ) awaits only a discussion of the area S a 
over which the integration is to be made. This discussion is important 
enough, however, to merit consideration in some detail and will be given 
in the following section. 


The Acoustic Plan Form 


Suppose that a line of sources is placed along the y axis and that 
the strength of these sources is zero for t < 0. At t = 0 they are 
"turned on" and, at the same time, start moving along the negative x x 
axis with the velocity V 0 . After time t* has passed, the source line 
has traveled a distance J^t as shown in the accompanying sketch (which 
is drawn for the case Mq>1). 


g 

Quotes are used around the word time since the dimension of t is 
actually length, not time. It is convenient, however, to refer to t 
as "time," and, since the actual value of time is simply t multiplied 
by the constant a Q , this should cause no confusion. 
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Suppose next that there are two sensing elements, or detectors, placed at 
the point P(x,y) located somewhere ahead of the y x axis; one of these 
detectors is responsive to light and the other to sound. Wow, the light 
detector at any given time will show the sources lying in a straight line 
Just as they would appear visually at each particular instant. The situa- 
tion is entirely different, however, for the sound detector. First it is 
necessary to understand the nature of a spherical sound wave. Such a 
wave travels outward from its origin at a velocity' a 0 , so that in the 
time t* it has traveled a distance t. Before the wave reaches a point, 
the point is completely unaware 4 of its existence and, further, after the 
wave has passed, the point remains subsequently unaware of its existence. 
Hence the only points disturbed by the wave are those momentarily on the 
spherical surface itBelf. (in this connection see reference 1, pp. 1 ff.) 


4 

Unaware is used in the sense that an instrument will record no change 
in any of the physical properties of the air at the point in question. 
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The sound detector, therefore, can only "hear" sources which are so 
located that their spherical sound waves are Just, at the given instant 
of time, reaching the detector. The locus of all the points which, at a 
time t ago, emitted sound waves that are Just now reaching the point 
p(x,y) is itself a sphere and for convenience this sphere will be 
referred to as an inverse sound wave. 5 The traces of these inverse sound 
waves in the z = 0 plane are drawn in the sketch as concentric circles 
about the point P(x,y). The intersection of an inverse sound wave of 
radius t— t with the line representing the position of the sources at a 
time t ’ — ' t* ago gives the position of the sources which are Just now 
signaling their presence to the sound detector at P. For example, when 
the source line started, it was lying along the yi axis. With reference 
to the present time t this was a 0 t ' removed. Hence the intersection 
of the circle about P of radius t with the yi axis fixes the two 
points A and A', the spherical sound waves of which are now reaching 
P. A continuation of this process yields, for the locus of all points 
from which waves emanated that are just now touching the point P, the 
part of the ellipse shown in the sketch. 

In the sense that the light detector, because of the very large 
velocity of light, is "seeing" a straight line of sources, the sound 
detector, because of the relatively slow velocity of sound, is "hearing" 
an elliptic line of sources. Extending this concept to include a sheet 
of sources distributed over the surface of the wing, one can refer to 
the outline of that part of the wing which generates disturbances which 
can be measured by the light detector as the plan form (i.e., the visual 
plan form), and to the outline of that part of the wing which affects 
the sound detector as the acoustic plan form. In a mathematical sense, 
the acoustic plan form is the area S over which the integration of 
equation ( 13 ) is to be made. 

The equation for the acoustic plan form can be formulated by means 
of the two equations 


(x— x x ) 2 + (y-yi) 2 + z 2 = (t--r) 2 (14) 


f(yi,Xi,T) = 0 


(15) 


where 

x,y,z coordinates of the point at which the induced effects are to be 
measured 

5 The inverse sound wave has for its analogue in steady lifting-surface 
theory the Mach forecone. In that theory a disturbance outside the 
Mach forecone cannot affect the values of the induced velocities at 
the point where they are being measured. Similarly, in the present 
study, a source located outside the inverse sound wave of radius t 
cannot affect the values of any measurement made at the point P. 
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variable points of the sources 
t "time" now 

t-i "time" ago 

Equation ( 14) is that of the inverse sound waves and equation (15) repre- 
sents the position of the visual plan form at a time t— r past. It is nec- 
essary to include the region behind the wing covered by the vortex wake 
as part of the visual plan form. In case the vorticity in the wake 
vanishes, as in the thickness problem, the wake may still be considered 
as part of the visual plan form, but the strength of the source-doublet 
distribution over that part of the acoustic plan form corresponding to 
the wake will vanish. If t 2 -z 2 <0, the acoustic plan form does not 
exist since, for such a case, no source on the *i>yi plane has had time 
to transmit its effect to the point x,y,z. On the other hand, if the 
circle (in the Xj^yi plane) given by the equation 

(x-x 1 ) 2 + (y-y 1 ) 2 =t 2 -^ 2 


lies entirely within the area occupied by the visual plan form of the wing 
at the beginning of the motion, the acoustic plan form is just this circle 
itself. For any other situation the acoustic plan form is formed in part, 
or in whole, by the curve found by eliminating t from equations (l4) 
and (15) and in part, or not at all, by an arc of the circle (x-Xi) 2 + 

( y— yi ) 2 = t 2 -z 2 . One of the principal advantages of the acoustic— plan- 
form concept arises in the study of problems involving source or doublet 
distributions having constant strength. In such cases, the retarded 
values of the potential and its gradient appearing in equation (13) are 
constant and can be. taken outside the integral signs. The problem is 
thereby reduced to the integration of a simple geometric variable over 
Sq. 


A few examples will serve to fix the idea of the acoustic plan form. 
Consider first a two-dimensional, unswept wing moving at a constant super- 
sonic speed in the negative x x direction. At time zero the leading edge 
of the wing was along the y x axis and now, at time equal to t', the 
wing has moved so that the leading edge coincides with the line Xi = — 
Mot. Choose three points that are now lying on the wing. Let one point 

have its x coordinate in the range c Q — ,/t 2 — z 2 > x > 2 (where c Q 

is the chord of the wing), the second in the interval J t^-z 2 >x > - 

J t 2 -z 2 , and the third is the range -,/t 2 -z 2 ^x ^ -Mot. Designating 
these points by P^Ps, and *3 (see sketch), it can be shown that their 
acoustic plan forms are, respectively, a complete circle, a part circle and 
part ellipse, and a complete ellipse. 
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The points P are at the centers of the circles and. at focal points of 
the ellipses. Since, moreover, the circular plan form about P x receives 
no signals from sources on the leading or trailing edge, conditions at Pi 
are consequently completely independent of the actual (visual) plan form 
of the wing. The elliptical plan form about P 3 , on the other hand, 
depends entirely on the shape of the leading edge; and finally the mixed 
plan form about P 2 is in certain regions (the circular portion) independ- 
ent of the leading edge, and in other regions (the elliptic portion) 
entirely dependent upon it. Since the wing is traveling at supersonic 
speeds, the trailing edge and vortex wake can have no effect on the meas- 
urements taken on the wing and , in the same way , a point ahead of the wing 
leading edge, P 4 in the sketch, is undisturbed. 

Next consider a wing moving at a constant subsonic speed in the nega- 
tive x x direction. As before, the leading edge was on the yi axis at 
t f = 0 and has traveled a distance -M 0 t. Choose now three points Pi, 

P 2 , and P on the wing and unaffected by the wing tips. The acoustic 
plan forms for these points are combinations of circles and hyperbolas as 
contrasted with the circle— ellipse combination in the supersonic case. 

Just as in the supersonic case, however, there is a certain region 
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represented by Pi in which the acoustic plan form is a complete circle 
and is independent of the visual shape of the wing (see sketch). 



*1 


1 

Point P 2 is surrounded by a plan form which is part hyperbolic and part 
circ ular , the point itself being the center of the circle and the focus 
of the hyperbola. Point P 3 is a limiting value of P 2 ; it lies on the 
leading edge of the wing and the hyperbolic sides of its plan form have 
degenerated into straight lines. Finally, P 4 lies ahead of the wing; 
its plan form is still a combination of a hyperbola and a circle, but P 4 
is now the focal point lying ahead of the hyperbolic branch used. 

The sketch was constructed so that the portion of the visual plan 
form behind the trailing edge had no effect on the potential at the var- 
ious points Pi, etc. If these points had been chosen at positions where 
the wake could signal its effect, one of two acoustic configurations would 
result. First, if the wing is symmetric about the z = 0 plane, no lift 
is developed and the vorticity in the wake is zero so that the visual plan 
form need not include the wake, but effectively ends- at the trailing edge. 
In this case, the leading edge of tho acoustic plan form is then deter- 
mined as before, while its modified trailing edge may be made up, in part. 
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of circular area formed by the primary wave and, in part, by an arc of the 
hyperbola formed by the (acoustic) intersection of the straight visual- 
trailing edge with the primary wave (such an arc being identical with the 
leading edge of the acoustic plan form but displaced backwards). On the 
other hand, if the wing has no thickness but is inclined to the free 
stream, it develops lift and the vorticity in the wake does not vanish; 
the acoustic plan form has a trailing edge made up entirely of an arc of 
the primary inverse sound wave. The space between this arc and the acous- 
tic trace of the visual trailing edge is covered by a sheet of doublets, 
the strength of which is determined by the vorticity distribution of the 
vortex wake . 

It is interesting to notice the conversion of terminology which 
arises in the analysis of unsteady lift problems. In the study of steady- 
state wings, it is customary (because of the nature of the governing 
partial differential equation) to speak of the subsonic problems as ellip- 
tic and the supersonic problems as hyperbolic. Yet the acoustic plan 
forms Just presented involved ellipses for the supersonic wing and hyper- 
bolas for the subsonic case. 

To complete the remark, it can be observed that when the velocity of 
the wing is sonic the steady-state partial differential equation becomes 
parabolic and, in this case, the acoustic plan form of a straight— edged 
wing also involves parabolas. When the leading edge is linear and normal 
to the stream direction, the eccentricity of the conic sections bounding 
the acoustic plan form for a point on the wing is equal to 1 /Mq. From 
this relation it is apparent that for Mq less than, equal to, and greater 
than 1 the sections are, respectively, hyperbolas, parabolas, and ellipses. 
As might be presumed, the value of the eccentricity satisfies simple sweep 
theory so that, for an infinitely long straight leading edge, the eccentric- 
ity of the acoustic plan form is 1/Mq cos A where A is the angle of 
sweepback. The principal axes of the conic sections are always normal 
and parallel to straight leading edges. 


Homogeneous Boundary— Value Problems 


Kirchhoff is solution to the wave equation can be applied to arbitrary 
wing plan forms undergoing arbitrary maneuvers. The boundary values for 
such general problems, however, usually lead to the development of double 
integral equations which are difficult to solve. As is usual in such 
cases, there are many special types of plan forms and maneuvers which lead 
to boundary— value problems that are simpler to analyze. An important 
class of these simplified problems is that arising from homogeneous bound- 
ary conditions. 

Let (p(x,y,z,t) be a solution to equation ( 3 ). I n certain special 
cases this can be written cp = (t) n <p Q (^r, ^r, ^) in which case cp is 



MCA TN 2256 


15 


called a homogeneous function of degree n. The number of variables 
affecting q> Q is only three as compared to the four which are necessary 
to determine q>. If, therefore, a partial differential equation can he 
set up for (p 0 , it will contain one less mathematical "dimension" than 
the equation for <p . Following this observation it is necessary to pro- 
ceed in two directions; one to find the partial differential equation for 
q> 0 , and the other to find the physical problem and consequent boundary 
values leading to a homogeneous flow field. The latter path will be 
first explored. 

First, consider an example 
of a homogeneous boundary— value 
problem. Suppose that a rectan- 
gular flat- plate at an angle of 
attack starts suddenly from rest 
and moves forward at a supersonic 
Mach number M 0 . At "time" ti 
the initial spherical wave gen- 
erated by the forward right-hand 
comer has traveled outward to a 
radius t x and, at "time" 2ti, 
to a radius 2ti. The sketch 
indicates the traces of these 
spheres in the z = 0 plane 
together with the original and 
present position of the wing 
leading edge. Let the points Pi 
and P 2 be located on the same 
rays through the origin of the 
circles and the wing comers. The 
problem is to find the pressures 
at Pi and P 2 . 

It is apparent that, if 
every dimension in the figure 
involving P 2 is divided by 2ti 
and every dimension in the figure 
involving P x is divided by tj., 
the two figures will be similar in every respect and point P x will coin- 
cide with point P 2 . Since the vertical velocity w 0 is constant over 
the plan form, a simple change in scale has made the boundary conditions 
for both problems identical. But this means that the solutions at Pi 
and P 2 are identical since the wave equation is invariant to change in 
scale. Hence, in regions of a rectangular wing unaffected by the waves 
from the trailing edge, the pressure can be written 



Ap 

q 



(16) 
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and the pressure is a homogeneous function of degree zero. A generaliza- 
tion of this example is contained in the following statement: 

(1) The pressure in any region affected by only two intersecting 
edges of a straight-sided flat plate traveling at a uniform 
subsonic or supersonic, speed is homogeneous and of degree 
zero (i.e., satisfies equation (l6)). 

Consider as another example the case of a flat rectangular wing 
traveling forward at a subsonic or supersonic speed and rolling about 
one edge, taken to be coincident with the x axis. The argument follows 
the same lines as before, and again a change in scale, proportional to 
the time, makes the geometry of the wing— wave combinations identical (in 
regions affected by only two intersecting edges) for different times. 

The boundary values over the wings will not be the same, however, unless 
the slope of the w 0 distribution is adjusted in each case. But w 0 
can be adjusted by' reducing it an amount proportional to the distance 
from the axis of rotation. The boundary-value problems are then similar 
for different values of time. Finally, therefore, the pressure can be 
written 


Ap 

q 


(x,y,z,t) 

\ 


= 7 


Ap 

q 



( 17 ) 


which is a homogeneous function of degree one. A generalization of this 
example is expressed as follows: 

(2) The pressure in any region affected by only two intersecting 
edges of a straight-eided flat plate traveling at a uniform 
subsonic or supersonic speed and rotating at a constant rate 
of pitch or roll is homogeneous and of degree one (i.e., 
satisfies equation ( 17 )). 

It should be noted that both (l) and (2) are equally true for the 
steady-state case when all transient effects have disappeared. In super- 
sonic wing theory they lead to conical and quasi-conical flows, respec- 
tively, while in the subsonic case they lead to flows about wings having 
infinite chordwise extent. In general, homogeneous flow occurs when the 
boundary conditions after a change in scale are proportional to their 
original values. 

Consider next' the modification of the basic partial differential 
equation (equation ( 3 )) under the assumption that the flow is homogeneous 
If the pressure is given by a function that is homogeneous and of degree 
zero, then, by equation (6), the velocity-potential function will be 
homogeneous and of degree one. If the notation 
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x y _ z _ 

— x-0 9 t - yoj t - zo 


<P (x,y,z,t) = t« (xq, y Q , z Q ) 


is used, then equation (3) becomes 


( 18 ) 


(l-x 0 2 )» x x + (l-y 0 2 )®y 0 y 0 + (l-^o 2 )*z 0 z 0 - 2x 0 y 0 $ x 0 y 0 


2x o z o ®x 0 z 0 2 *o z o $ y 0 z 0 0 


and a linear partial differential equation with three independent vari- 
ables is therefore obtained. 


In the general theory of 
partial differential equations of 
second order, the character of an 
equation is determined from the 
geometric nature of a related 
quadric surface. The character of 
equation (19) can be shown from 
such considerations to depend on 
the sign of the expression 1 — Xq 2 — 
y Q 2 -z 0 2 . It is immediately 
apparent, however, that within 
the unit sphere in the x^y^z,-, 
space the sign of 1 — Xq 2 — yo 2 ~Zo 2 
is everywhere positive and outside 
the sign is everywhere negative. 

It follows that outside the unit 
sphere equation (19) is hyperbolic 
and inside the unit sphere it is 
elliptic . 

The character of equation (19) 
is of particular interest since the 
difficulties inherent in the deter- 
mination of the solutions can be 
estimated without actually obtaining 
the solutions. For example, con- 
sider the two configurations shown 
in the sketch. These wings started 
moving at t = 0 with the foremost 



portion of their leading edges on 

the y 0 axis and have by now traveled forward at a supersonic speed to 
attain the positions represented by the sketch, the unit circle being in 


each case the trace of the primary wave from the vertex on the z = 0 
plane. Outside the unit sphere, the governing equation is hyperbolic and 
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the behavior of the flow is similar to that in steady-state supersonic- 
wing problems; closed solutions, therefore, can be expected in many cases. 
Inside the unit sphere, on the other hand, the character of equation (19) 
is elliptic. The canonical fonn of an elliptic partial differential equa- 
tion in three dimensions is Laplace’s equation which is the governing 
equation in steady-state subsonic-wing theory. Judging from the com- 
plexity of subsonic, three-dimensional wing problems, it can be concluded 
that a closed solution within the unit sphere will be very difficult to 
find. 


It is instructive to notice that this entire development has a direct 
analogue in the study of three-dimensional, steady state, supersonic wings. 
In that case the original equation is the three-dimensional wave equation 


'Pxx-'Pyy— ^7= 0 (20) 

By considering the velocity potential to he homogeneous and of degree one, 
Busemann in reference 7 was able to introduce the transformations 



<p(x,y,z) = x4>(y 0 ,z 0 ) 
and transform equation (20) to the form 


(l-7o 2 ) ®y 0 y 0 + (!-Zo 2 ) $z 0 z 0 -2y 0 z o ®y 0 z 0 = 0 ( 21 ) 


which is the two-dimensional form of equation ( 19) . Flows governed by 
equation (21) have become known as conical flows. A study of equation (21) 
shows it to be elliptic inside and hyperbolic outside the unit circle. 

In this case, however, the equation has only two independent variables so 
that once the equation has been transformed to the two-dimensional form 
of Laplace’s equation 6 solutions^ are not difficult to find. 


0 

The Tschapligin transformation transforms equation (21) into Laplace’s 
equation. 
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Boundary— Value Problems Involving Noninteracting Surfaces 


The simplification of the four-dimensional theory brought out by 
the introduction of homogeneous flow was more apparent than real since 
the resulting partial differential equation, although containing one 
less dimension, was unwieldy. Another class of wing problem which is 
simplified both in theory and in practice by reasoning from physical 
knowledge of the flow behavior is that in which the wing has a super- 
sonic edge (i.e., an edge which is traveling with a supersonic normal 
component of velocity). 

When the acoustic plan form is affected only by a supersonic edge, 
it is obvious that the flow on the upper surface of the wing is independ- 
ent of that on the lower surface. Hence the solution to such problems 
can always be written in terms only of sources as follows: 


""-ST. 


1 d(p 

r o LszT. 


I d a=0 

(b) 


a^(t). 


where dcp/dzi = w u (xi,yi) = 

V 0 X u (xi,yi), Xu being the T V///// 7 '' '''/////// 

local slope of the surface in c v/y/a=<x(t) '/////// 

the direction of V Q . Since _i_ Y/// /// / /////// / /X 

equation (22) is equally valid ( a ) 

for symmetrical nonlifting ^ W / ////// ( 

surfaces and lifting plates, c y a =a(t)/A a=0 \ 

its value and simplicity is Jl- {// / / / / / A « 

evident. (b) 

If the wing plan form is "T" a =0 

further specialized by having ^ , K/s / / s / / /\ 

not only supersonic leading ( c ) - - - - - - 

edges, but also having a ___ . f— q.=<xv, ; _ _ 

straight trailing edge perpen- | X 

dicular to the direction of ^ ^ 

motion, additional simplifica— — A 

tions can be used. 7 Consider, ' ' r a=a ' ' 

for example, the two-dimensional 

wing (a) in the sketch. Let — ^ yT 7 ] 

this wing have an angle of c x Va a=0 x 

attack <x(t) which varies with ^ / 1X1 \ 

time in an arbitrary manner. 

There results from such an 
angle— of-attack variation a 
certain lift which also varies 

with time. Hence, if L* represents the total lift on an airfoil of 

7 The following method simply extends, to include the effects of unsteady 
motion, a theorem given by Lagerstrom and Van Dyke. (See reference 2.) 
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very high aspect ratio and * 
slope, then 8 a 


represents the section lift-curve 


— = c c 7 * (c,t)(span) 
qa a 


Next, it is clear ty reason of symmetry that the total lift on wings (b) 
and (c) in the sketch are equal. Then, since the analysis is based on a 
linear partial differential equation, by superposition principles the 
total lift on wing (b) or (c) is Just half of that on wing (a). In 
another sense, the lift coefficient for the whole wing based on the 
deflected area is the same for all three cases. A suitable superposition 
of wings (a), (b), and (c) will give wing (d), which then has the same 
lift coefficient based on deflected area. Finally, because of the super- 
sonic stream, wing (e) can be obtained from (d), hence it also has the 
lift coefficient common to the other wings. It is of course necessary 
that the variation of a with time be the same in each case. 

The preceding process can be extended one step farther to the 
development of the lift due to a single deflected element. By consider- 
ing the sketch it can be seen that 



L^_ 

qa 


c 


c 2 a '(c,t) 


(c-£c)c 7 ' (c— £c,t) 

AS 

Ac 


where 

AS area of the deflected element 
c distance from AS to trailing edge 

By the usual limiting process the latter equation becomes 


L»_ _ _S_ 
qa Sc 



(c,t) 


dS 


e The primes on quantities indicate that two-dimensional values are taken, 
or, that a high-aspect-ratio wing is considered and tip effects are 
neglected, • - 
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Finally, if a wing is composed of a distribution of these elements, then, 
for an axial system centered at the apex of the wing leading edge, there 
results 

C L =— | f J a(x,y) ^ (c D -x) c la » (c Q -x,t) J dydx (23) 
s 

where Co is the maximum chord (aee sketch ) . In the development of equa- 
tion (23) each element is assumed 
to have the same variation of 
motion with time. 

Notice that when the wing is 
a flat plate flying at a steady 
speed so that all transient effects 
have disappeared, cj 1 is independ- 
ent of t and of the chord length, 
being, in fact, equal to k/$. Then 
equation (23) becomes 

° L ' ikj. / a<I ’ y)toljr = T? 

where a is the average angle of attack of the wing; This result has 
already heen obtained in reference 2. When <x(x,y) is independent of 
x (as for a flat wing sinking or rolling), equation (23) becomes 

b 

C L = c /1W) c C7 ' (c,t) dy (24) 

o Jo Q> 

2 

vhere c is the local chord which is, in general, a function of y. 

Squat! on (24) simply indicates that longitudinal— strip theory is exact 
Tor such wings. 

Finally, notice that the calculation of the unsteady lift on three— 
limensional wings with supersonic leading edges and straight trailing 
sdgos perpendicular to the free-stream direction has been reduced to an 
integration involving the relatively simple results for a two-dimensional 
■ring undergoing the same unsteady motion. For example, the lift on a 
flat triangular wing rising and sinking with a harmonic motion can be 
computed from a single integration of the results presented in. refer— 
snce 3. Such a calculation can be carried out numerically quite rapidly. 

Two-Dimensional Boundary-Value Problems 

The simplification brought about when the flow is independent of 
•ne dimension is again obvious. In such cases, the three-dimensional wave 
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equation (3) reduces immediately to the two-dimensional wave equation. 
Typical examples of this type of problem can be constructed by consider- 
ing flat plates which start suddenly at t = 0 and travel thereafter at 
constant supersonic velocities. Two examples, one a corner of a rectan- 
gular wing and the other a triangular flat-plate wing are shown, in the 
sketch. After time t = 0, the edges of the wings send out cylindrical 
waves and the outer boundaries of these waves at time t are shown as 
dashed lines parallel to the edges in question. 



Since points in regions 1 and 2 are affected only by a single edge, the 
wave phenomena in these regions are cylindrical, and the physical quan- 
tities are in both cases independent of distance parallel to the edge 
which acts as their generator. Hence, the flow field in these regions 
may be regarded as two dimensional. (Region 3* incidentally, is independ- 
ent of distance in both the x and the y directions and is, therefore, 
one -dimensional. ) 

Solutions to the two-dimensional unsteady problems are sometimes 
especially easy to find because of the analogy they have with three 
dimensional, steady-state, lifting-surface problems. (See, e.g., refer- 
ence 8,) For example, consider an infinitely long unyawed wing which 
starts from rest and travels forward at a velocity V which may or may 
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not be a function of time. The trace of 
this wing in the x,t plane is like that 
shown in the .sketch, (in the sketch shown, 
the wing velocity is always less than the . 
speed of sound and varying.) The boundary 
conditions are that <p z is specified over 
the shaded area and the loading A<p^ is 
zero everywhere except within the shaded 
area. But if x is replaced by y and 
t by x, these boundary conditions are 
exactly the same as those for a plate of 
known camber and angle of attack, with a 
plan form as indicated by the shaded area, 
placed in a free stream directed along the 
positive ' x axis at a Mach number equal 
to y2. The solution for the one problem 
may be used, therefore, as a solution to 
the other with only a change in notation. 


Traces of 


characteristic cones 



Boundary Conditions for Very Slender Wings 


When the wing plan form is slender in the sense that its length in 
the streamwise direction is large compared to its span, an estimation of 
the loading on it can be 
obtained by neglecting in 
the partial differential equa- 
tion the gradient of the induced 
velocity component in the stream 
direction. Thus, if the wing 
is moving in the negative x 
direction, equation (3) reduces 
to 


< * > yy + c * ) zz ^tt 


( 25 ) 


which is again the wave equa- 
tion but in two space dimensions. 
Since equation (25) is independent 
of x, study can be made independ- 
ently of the boundary conditions 
on each span wise strip. This is 
an extension of steady-etate 
slender-wing theory. (See, e.g., 
reference 9.) The sketch shows a 
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typical section in the yt plane. If the wing is a flat plate starting 
from rest and traveling forward at a uniform speed, the boundary condi- 
tion is that <P Z should be constant over the shaded area in the sketch 
and A<p+. should be zero everywhere except across this area. The solution 
to this problem leads directly to the solution for the thin triangular 
wing. 9 


Two-Dimensional Unsteady Incompressible Flow 


The analogy between two-dimensional unsteady and three-dimensional 
steady flow provides an interesting viewpoint for two-dimensional unsteady 
problems in incompressible flow. In this case the boundary conditions 
are independent of the lateral coordinate and equation ( 1 ) becomes 

1 

Txx+'Pzz = —z 'Pft* 
a o 

Since the speed of pressure propagation a 0 is assumed infinite, the 
basic equation thus becomes 


<Pxx-*Pzz = 0 



where, however, time still appears 
in the boundary conditions and in 
the expression 

Ap _ Jl_ l£_ 

q " V 0 2 at* 

which is equivalent to equation (6). 
The technique involved in the solu- 
tion of such problems is analogous 
to that in which one term of the 
three-dimensional steady— state 
equation vanishes by virtue of the 
fact that free-etream Mach number 
approaches one. 

If the two-dimensional wing 
starts from rest and travels for- 
ward at a speed V, the trace of 
the wing is as shown in the sketch. 


The essential difference between this problem and the more general case 
of two-dimensional compressible flow lies in the fact tha t the traces of 
9 It is obvious that the analogous problem in steady-state wing theory is 
that of a low-aspect-ratio, rectangular, flat plate in a free stream 
having a Mach number equal to ^2. . 




NACA TN 2256 


25 


the characteristic cones are normal to the t’ axis. The boundary 
conditions axe therefore satisfied along lateral strips and, in lifting- 
surface terminology, the. analysis corresponds to slender-wing theory.. 
These latter methods are well established and in reference 9, for example, 
the manner in which the Kutta condition is imposed is discussed in some 
detail. The trailing vortex sheet for the lifting wing has the same 
distribution of vorticity that exists behind the unsteady two-dimensional 
airfoil and the rolling up of the vortex sheet can be studied from either 
standpoint . 


APPLICATIONS 


Starting Lift of a Wing 


One of the simplest and yet most general results which can be derived 
on the basis of the present theory is the initial value of pressure on a 
wing surface starting impulsively from rest with a velocity V 0 . The 
discussion will be made for a wing without thickness although the results 
will be seen to apply to the thickness case as well. : 


Consider a surface with a 
plan from as indicated in the 
sketch. The acoustic plan form 
of a point P(x,y) on the sur- 
face is a small circle of radius 
t. Since no point on the wing 
outside this circle can influence 
the pressure at P, the upper 
surface is independent of the 
lower surface, except for a band 
of width t around the edge of 
the wing. It is, therefore, 
evident that the boundary— value 
problem to be solved has been 
treated in the section Boundary- 
Value Problems Involving Noninter— 
acting Surfaces. The solution 
follows directly from equation (22) 
and can be written 


1 nr 1 dq> 

2 it J J a r o Uzi 

b a 


dS 
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Using a polar coordinate system defined by 


x — Xx = r cos 6 
y - yi = r sin e 
dxi dy x = rdrdQ 


there results 10 


- 


w u (x,y) 


2 it ot 


2 ji 


r- x 


dr = — tw u (x,y) 


so that 


^) t=0 ” + vjloC^) 


t =0 


4w u (x,y) 

V0M0 


If a = - w u (x,y)/V 0 is the local slope of the wing, the egression for 
load coefficient becomes 


( ££) = ( 26 ) 

4=0 


The starting value of lift coefficient can, therefore, be written 

cl - jg ««) 

where a is the average angle of attack of the surface defined by 

S.iy'jTadxdj 

S being the area of the wing plan form. 


10 The mean value theorem gives 



where 4 and T| lie somewhere in 8&. Hence, as t approaches 0, 
4 and t| approach x and y, respectively. 
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Supersonic Steady— State Lift 


In order to illustrate some of the basic concepts, consider next 
an example for which the result is well known. The problem will be to 
derive, by the methods outlined in this report, the steady— state loading 
on a two-dimensional flat plate traveling at a constant supersonic speed. 

Since the upper and lower surfaces are noninteracting, the solution 
stems from equation (22). In the plane of the wing this equation becomes 



w u dx x dy x 
(x-x x ) 2 +(y-yi) 2 


( 28 ) 


Leading edge now 



The accompanying sketch shows the positions of the wing in the xy 
plane. The wing has constant speed for t>0. The point P(x,y) is 
chosen on the wing and ahead of the wave which started at time zero; 
therefore, P(x,y) lies in the region which has attained its steady-state 
value. Further, the value of is constant over the acoustic plan 

form in such a region. This constancy reduces the problem to one of 
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integrating [(x-xx) 2 +(y-yx) 2 ] x / 2 over the ellipse representing the 
acoustic plan form. 

The equation for Sg can .be determined from equations (14) and 
( 15 ). In this case, equation ( 15 ) becomes simply 

xi = - MqT 

Elimina te t between this equation and equation (l4) and there results 
in the z = 0 plane 

(x-xx) 2 + (y-yi) 2 = (t + ( 29 ) 

That equation (29) is the equation of an ellipse with one focal point 
at x,y can be readily verified. It is more convenient, however, to 
change to a polar coordinate system 

x— x x = 

y-yi = 

dyi dxx = 

Then equation (29) becomes 

Mor = Mot 


Mq + cos 9 


with origin at P. Hence set 

r cos 6 
r sin 6 
rdrd8 

+ x— r cos 9 


+ Mot 


( on \ 


and, therefore, equation ( 28 ) becomes simply 

x + Mot 

, . _ jjjl 4 008 8 dr 




de 

+ cos 9 


The integral is not difficult to evaluate so 


<P 
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and, finally, by equation (6) 

Ap _ 4 f _ v uMq ^ _ 4a ^ 

« = v oMo v 

which is the familiar Ackeret value for the loading on a two-dimensional 
flat plate. The lift coefficient, of course, follows immediately as 


C 


L = 


4a 

A> 2 -i 


(32) 


Stability Derivative 


Consider the case of a high-aspect-ratio (so that the tip. effects 
can be neglected) rectangular wing yawing at a uniform rate r radians 
per second and traveling forward along its flight path at a supersonic 
velocity V 0 . The equation of the leading edge of the wing, using the 
coordinates indicated by the sketch, is given by the relation 


x 





and since = rt* this 
becomes 

tan rt‘ (33) 

The boundary conditions 
(when the stability derivative 
is to be referred to the body 
axis) are that w u , the ver- 
tical induced velocity in the 
region occupied at any given 
time by the plan form, shall be 
constant. Further, the steady- 
state condition is determined 
when the intersection of the 
leading— edge trace with the 
inverse sound wave forms a 



Wing 



closed curve. Honce, the 

acoustic plan form is given by eliminating t from the equations 
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(x-Xi ) 2 + (y-yi ) 2 = (t-T ) 2 
xi = (yi - tan r-r' 


( 3*0 


It is usual to assume r to be small and. to retain only the first 
power of r in a series expansion. Therefore, since it can be shown 
that the equation of the acoustic plan form is given correctly to the 
first power in r by assuming rt' = tan.rt', equation (3*0 becomes 


R = t — T 

x — R cos 9 = ( y — R sin 6 


=(7- 



(35) 


where, as in the preceding example, the polar coordinates 


x—Xx = R cos 9 
y-yi = R sin 9 
dx x dyi = RdLRdQ 

have been introduced together with the notation 



solving equations (35) f° r there results 

-A 0 - U B 0 r* sin ~9 

2r* sin 9 


( 36 ) 


( 37 ) 


where 


A q = cos 9 - r*y + Mq - r*t sin 9 
B 0 = r*ty— x-Mot 

and expanding in powers of r* gives, neglecting powers higher than the 
first. 




r* sin 9 


or 


R = 


x+MQt 
Mq+ cos 9 


- r 1 


ty 


(x+Mot)(y+t sin 9) (x+Mot)' 


LMo + COS 0 (Mq+ cos 0)‘ 


(Mo+ cos 9) 3 


sin 9 
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The solution for the velocity potential is given, just as in the 
preceding example, by the equation 

^ 2JT 




and since 


r 21t sin 0 _ r 2n sin 0 

l Q (Mq. + cos 0) 3 J 0 (Mq + cos 0) 2 


p 

J ~ Mn + 


Mq + cos 0 VMq 2 -1 


/ 2Jt 

(mT+ 


d 0 _ 2 nMp 

(Mq + cos 0) 2 “ (Mq 2 — l) 3 / 2 


then finally 


= _ vpsC" 0 * + 1 * r,y 


The loading can be obtained from equation ( 38 ) by using equation ( 6 ) 


hence 


Ap _ 4a, 

<J “ V Mq 2 — 1 + (Mq 2 -!) 


\ w 


This result has "been obtained previously in reference 10 by a different 
method. 


Moving Two-Dimensional Vortex 


Moving 

vortex 


A further interesting application of the Starting 

methods presented here is the development of ^ vortex 

the moving rectilinear vortex, where a vortex ^ / 

is represented by a line of constant circula— T* w “ 

tion. The moving vortex is associated with a * *"j 

lifting surface 11 which starts impulsively at 
t = 0 and travels at a uniform velocity away 
from a fixed axial system. Since a vortex 

cannot end in a fluid (Helmholtz's fundamental result), the moving vortex 
must he adjoined by two trailing vortices and they, in turn, by a final 
starting vortex (see sketch) which remains in the vicinity of the position 

11 Associated in the sense that a distribution of these vortices will 
satisfy the boundary conditions for the lifting surface. The moving 
vortex is often 
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at which the motion commenced. For the present, the assumption is made 
that the trailing vortices are far enough away so that they do hot affect 
the flow arodnd the center of the moving and starting vortices ( i • e • , the 
flow is two dimensional); and second, that the starting vortex remains 
fixed on the line along which the moving vortex first appeared. With 
these assumptions, the boundary conditions can be determined from the 
definition that the vortex system described is equivalent to a distribu- 
tion of doublets of constant strength over the region of the z = 0 
plane bounded by the moving vortex, the stationary vortex, and the vor 
tices joining them at their extremities. Such a sheet of doublets gives 
rise to a constant potential difference, A 9 , across the sheet, while the 
normal derivative is continuous so that = 0. The solution is 


then given by equation (13). Thus 



9 = 9i + 9 2 


(40) 

where 



<Pi ■ 

| dS 

(41) 

and 



ifSss') 

knJJc r 0 \ dt J 

^ a 

dr 0 

— - cLS 

dz 

(42) 




The evaluation of 9i imposes 
no difficulties except in integra- 
tion but the calculation of 9 2 
requires some special consideration. 
Since A 9 is a constant, dA9/dt 
is zero everywhere except along 
the line of the moving vortex and 
there it is infinite. In order to 
evaluate such an indeterminate 
form, assume that the potential 
difference rises linearly to its 
value A9 in the time interval e 
(see sketch). Then the value of 
dA9/dt is Ap/€ in the narrow 
strip between the two ellipses 
formed by the intersection of the 
inverse sound wave 

(x-Xi) 2 + (y-yi) 2 + z 2 = (t-r) 2 
with the lines 

= “ MqT 

Xi = - M 0 (t- e) 
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When the equations for these two ellipses are formulated and placed in 
equation (42) there results an equation which can be represented by the 
expression 


where 


«P 2 = ^ F(t)-F(t— e) 


b a 


dS 


and So is the same acoustic plan form used in equation (4l). But since 


lim 
e —>o 


F( t)-F(t— e) = 

e St 


the equation for 9 2 becomes 


4« St J J c r o Sz 


The calculation of 9i and <J > 2 over the acoustic plan form gJven 
by the area within the curve 


(x-xj .) 2 + (y-yj 2 + z 2 = ^t + ^ 


(43) 


or between this curve and the yi 
axis gives for the regions defined 
by the sketch 

Region (l) 


2 _ ,,2 , _2 


q> = £2 
Y 2« 


^ (x+Mot) yt 2 -x 2 -z 2 Moving 

arc tan : vortex 


z(MoX+t) 


arc tan Jt 2 —: 

Z t / 


2 2 
•X -z 


(44) 



X 

Starting 

vortex 


Subsonic 


Region (2) 


<p = 

Y 2 rt 


l/l*l 


arc tan 


Region ( 3 ) 
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It is to be noted that equation (44) gives the potential everywhere 
within the cylinder in which disturbances can be affected in subsonic 
flow as well as over a large portion of the disturbed region in super- 
sonic flow. 

It is apparent that the results given in equations (44), (45), and 
(46) can also be obtained independently by the methods indicated in the 
section entitled Two-Dimensional Boundary— Value Problems. From this 
viewpoint, the problem of the unsteady two-dimensional vortex pair has 
for its analogue in lifting-surface theory the case of a triangular wing 
yawed so that one side is parallel to the free-etream direction and the 
other side is either a subsonic or a supersonic edge, depending on the 
case being considered. Since A<p is specified a constant over the 
entire plan form, the solution is relatively easy. In reference 11., for 
example, the solution for the subsonic case is given in terms of the 
variables used in steady-state analysis. 

The foregoing results for f^><l reduce to more familiar expres- 
sions when their asymptotic expressions for large values of time are 
determined. A study of the two terms in equation (44) shows that as t 
becomes large the expression becomes 

Acp / n z . x \ 

f * 27 {2 in - arc z) 

This result provides the velocity potential for a stationary vortex in 
an incompressible medium, the jump in the value of the potential occurr 
j n g along the plane of the vortex motion (see sketch). 

If the axes are fixed in 
the moving vortex, new coor- 
dinates I, £,t need to be 
introduced as defined by the 
equations 

| = x + Mot 

t • z 

T = t 

As t becomes large, equa- 
tion (44) reduces to the form 

cp > — ( — — ,4i + arc tan r—====\ 

2 A 2 |£j £ -/L-Mq 2 / 

This result agrees with that of Glauert (reference. 12) for a vortex 
obtained from the steady-state linearized theory of compressible flow, 
so that in the vicinity of the moving vortex the induced velocity field 
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is identical to the perturbation field of a fixed vortex in a subsonic 
free stream. 


It is interesting to inspect the vertical velocity induced by the 
vortex pair in the z = 0 plane. From equations (44) and (45) it 
follows that for x 2 <t 2 the vertical velocity w u for both subsonic 
and supersonic speeds is given by the relation 


_ f" ( 1 ] 

Wu " L Vx+Mot x) J t 2_ x 2 J 

mq [ 

2« x{ x+Mot ) ^ J 


(47) 


The sketch shows the 
variation of w u along the 
x axis for a moving vortex 
traveling at both subsonic 
and supersonic speeds. For 
both cases in the vicinity 
of the starting vortex w u 
behaves as it does near a 
two-dimensional, steady-state, 
fixed vortex, having large 
values of downwash on one 
side and large values of 
upwash on the other. The 
velocity distribution around 
the moving vortex is quite 
different, however, depend- 
ing on whether its Mach 
number is greater or less 
than one. In the subsonic 
case the distribution is 
similar to that for the 
starting vortex, but in the 
supersonic case the velocity 
is just a pulse existing at 
the position of the vortex. 
Again in both cases the 
velocity falls to zero at 
the primary wave front (i.e., 
at x = ± t in the sketch). 



A final remark can be made concerning equation (44) in connection 
with incompressible theory. In this case, the value of a Q becomes 
infinitely large and, if equation (44) is written in the form 
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m A9 
* = 


arc tan 


( x+Vpt [ ) ya Q 2 t ' g-x 2 -^ _ 


arc tan 


x ya 0 2 t 1 g-xg-z 2 

za 0 t' 


^ + a o t0 

a o 0 y 

it is apparent that, as a 0 ->oo the velocity potential hecomes 


-g (• 


X-t-V 0 t* 

arc tan arc tan 

z 


f) 


(48) 


Equation (48) gives the velocity potential for a moving and a fixed 
vortex pair in an incompressible fluid. 


Moving Three-Dimensional Vortex 


In three-dimensional subsonic 12 unsteady lift theory an important 
role is played by the loop vortex composed of a moving bound vortex of 
fixed span, a stationary starting vortex of opposing strength, and two 

trailing vortices connecting the 
end-points of the other two and 
completing the vortex loop. The 
boundary condition for such a 
configuration requires that the 
jump in potential A9 is a con- 
stant over the interior of the 
loop. The calculation of the 
induced field can be resolved 
essentially into the determina- 
tion of the velocity potential 
for two semi-infinite vortices 
and a trailing vortex (see 
sketch) since, by superposition 
methods, two of these latter 
configurations can be combined 
with the two-dimensional figure 
of the last section to give the 
required loop. The present 
example is concerned with the calculation of the velocity potential in 
the tip region of the semi— inf ini te case. 

As in the two-dimensional case just discussed the potential can be 
divided into two parts 9i and 9 2 such that 


12 The extension to include the three-dimensional supersonic loop is 
obvious. 
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9 = 9j. + 9 2 


m , A9z p p ^i d yi 

1 " 4* JJ q [ ( x— Xi ) 2 +( y— yi ) 2 +z 2 ] 3 / 2 

b a 

m a = + £21. A r f teiAZll 

I+Jt dt j Jq [ (x-x 1 ) 2 +(y-y 1 ) 2 +z 2 ] 

d 

where S a is the acoustic plan form* being the area hounded by the 
lines 

(x-*i ) 2 + (y-yi ) 2 + z 2 = ^t + 

yi = 0 

ij = 0 


It is convenient to use the notation 

k = z+Mot 
I = X-X! 

ti = y-yi 

Then the potential can he written 


cp = 


zAcp p y 

kn Ayt 2 -x 2 -z 2 








2 + p 2 z 2 + p 2 T) 2 ) 


di 

U 2 + ^ + z 2 )3/2 


zAcp £ p y 

w st J r i 2 2 

— v t — x — z 



(k-M 0 yk 2 +p 2 z 2 +p 2 ri 2 ) 


dj 

i^+Tp+z 2 
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which becomes, after one integration, 


r y ' 

ffi = SM / 

kn J j 


2 2 
■X -Z 


-Mp X 

Lk 2 +Z 2 +11 s (ti 2 +z 2 )7x 2 +Z ^ 


2 +p 2 z 2 +0 2 


?/’ [ 

- /t 2 -: 0 0 L 


Mo‘ 


: +Z 2 k 2 +TJ 2 +Z 2 


■)] 


dq + 


Mo 




k 2 +z 2 +T) 2 (k 2 +z 2 +q 2 ) V^ 2+ P 2 z 2+ P 2t 1 2 


■X 2 — z 2 


dT) 


Combining and integrating again there is obtained, for Mq<1, 


cp = r 
4rt |_ 


( x+Mot ) J t 2 — x 2 -z 2 Xa/ t 2 — : 

arc tan - ii — ; arc tan 

z(MoX+t) tz 


2 2 


arc tan 


y(x+M 0 t) 

z ( x+M 0 t ) 2 +p 2 ( y 2 +z 2 ) 


— arc tan 


a 


xy 

x 2 +y 2 +z 2 


(49) 


Moving vortex 



which ia the potential within 
the primary spherical wave 
shown in the sketch. Notice 
that across the shaded area the 
jump in potential given by 
equation (49) is equal to A<P 
and elsewhere the potential is 
continuous. Along the portion 
of the sphere t 2 = y 2 +z 2 +x 2 
for which y> 0 , the potential 
is equal to that derived for 
the two-dimensional vortex; and 
on the rest of the sphere, 

(p equals 0 . 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., October 12, 1950 . 
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